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Abstract. We study, by numerical simulations, the en- 
trainment process in a supersonic, radiative jet flow, dur- 
ing the evolution of Kelvin-Helmholtz instabilities, in the 
context of the the acceleration problem of molecular bipo- 
lar outflows, observed in Giant Molecular Clouds. Our re- 
sults show that a large fraction of the initial jet momentum 
can be transferred to the ambient medium by this process. 
We therefore analyze in detail the instability evolution 
and compare some of the main observational properties 
of molecular outflows with those of the entrained material 
that we observe in our simulations. In particular, we find 
a good agreement for the mass vs. velocity distribution 
and for the outflow coUimation structure, especially when 
a light jet is moving into a denser ambient medium. This is 
probably the case for (obscured) optical jets driving pow- 
erful molecular outflows in the denser environment of the 
inner regions of molecular clouds. 

Key words: Hydrodynamics - instabilities - stars: early- 
type - ISM: jets and outflows 



1. Introduction 

Young Stellar Objects (YSO) are accompanied by different 
kinds of outflow phenomena, whose most striking man- 
ifestations are optical jets and molecular outflows; how- 
ever the mutual relationships existing between such dif- 
ferent kinds of outflows still remain to be clarified. In the 
past, the prevailing idea was that fast atomic winds, with 
low collimation, were the primary accelerating agents of 
molecular outflows (Lada 1985; Lizano et al. 1988), while 
optical jets were considered as parallel phenomena with no 
dynamical significance, because their momentum trans- 
port rate was estimated to be much lower than that of 
molecular outflows. These models, however, failed to ex- 
plain, in a satisfactory way, some important outflow prop- 
erties, while, on the opposite, models based on highly col- 
limated flows appeared to be able to account for them 
(Masson & Chernin 1992; Stabler 1994). Additional prob- 
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lems to the wind models have been posed by the discovery 
of molecular outflow components with a higher degree of 
collimation and higher velocities with respect to the stan- 
dard component, such as EHV components and molecu- 
lar bullets (Bachiller 1996). Moreover, the estimates of the 
momentum carried by optical jets have been reconsidered: 
more accurate comparisons of the observed emission with 
shock models has allowed to show that jets have a low 
degree of ionization (Raga, Binette & Canto 1990) and 
therefore their density is much higher than what was pre- 
viously thought and, secondly, the time scales of the jet 
phenomenon appear also to be much longer than what 
was previously thought (Parker, Padman & Scott 1991; 
Bally & Devine 1994); therefore the total momentum that 
optical jets can deposit in the ambient medium is much 
larger than what was previously estimated and the prob- 
lem of momentum deficit seems to be overcome. As a con- 
sequence of all these considerations, there is a growing con- 
sensus that there must exist a causal connection between 
optical jets and molecular outflows (Bachiller 1996), how- 
ever the physical mechanism through which jets drive the 
outflows have not yet been clarified. The proposed mod- 
els fall essentially into two classes: in one case the ambient 
material is entrained in a turbulent mixing layer along the 
length of the jet (sometimes this is referred to as steady 
state entrainment), while, in the second case, the entrain- 
ment is performed by a working surface, either the leading 
bow shock or internal working surfaces due to variability 
in the jet emission (sometimes this is referred to as prompt 
entrainment). 

While there are several studies concerning the second 
mechanism (see Cabrit et al. 1997 and references therein), 
for the first one there are only a few attempts to model 
the physical processes involved (Stabler 1994, Taylor & 
Raga 1995). Both mechanisms should be at work: work- 
ing surfaces and bow shocks are surely present in jets and 
the development of a turbulent mixing layer between jet 
and ambient material seems unavoidable as the result of 
the evolution of Kelvin-Helmholtz instabilities. Although 
in several situations the prompt entrainment mechanism 
seems to fit very well the observational properties (see e.g. 
Davis et al. 1997), a thorough study of the steady state 



M. Micono et al.: On the matter entrainment by stellar jets and the acceleration of molecular outflows 



entrainment is also needed in order to assess its impor- 
tance. In any case, more work is needed in the theoret- 
ical analysis of both mechanisms since they both seem 
to fail in explaining the large size of molecular outflows 
as compared to that of optical jets. In this paper we ad- 
dress the problem of steady state entrainment; its study is 
obviously complicated by the fact that it involves turbu- 
lence and one possible approach, which we will follow, is 
the use of direct numerical simulations. We have studied 
the evolution of Kelvin-Helmholtz instabilities in several 
different conditions both in two and in three dimensions 
(Bodo et al. 1994, 1995, 1998, Rossi ct al. 1997, Micono 
et al. 1998, 2000) and in this paper we analyze the results 
of three-dimensional simulations in a radiative situation. 
The study of a fully three-dimensional situation is essen- 
tial for analyzing the entrainment properties, since it is 
well known that the properties of turbulence are very dif- 
ferent in 2D and in 3D (for a comparison of 2D and 3D 
results on the evolution of Kelvin- Hclmholtz instabilities 
see e.g. Bodo et al. 1998). Following the above consider- 
ations, the focus of our analysis will be on the properties 
of the entrainment of the ambient material, while a more 
general discussion on the evolution of the instability is 
presented in Micono et al. (2000) (hereinafter referred as 
Paper I). 

The plan of the paper is the following: after a summary 
of the properties of molecular outflows, given in Sect. 2, 
we present our model in Sect. 3, and our results in Sect. 4. 
Finally, a summary and the conclusions are given in Sect. 
5. 



2. Properties of molecular outflows 

Molecular outflows consist of a pair of oppositely-directed 
and poorly collimated lobes, symmetrically located about 
an embedded YSO, detected in broad mm- wave emission 
lines, especially of CO. Their maximum sizes range from 
0.04 to 4 pc, the dynamical ages from 10^ to 2 x 10^ yr, II2 
densities are ~ 10^ cm~'^, their total momentum is in the 
range 0.1 - 1000 M© km s^^ (Bachillcr 1996, Fukuy et al. 
1993). In general, they are poorly collimated, with ratios 
of major to minor axis of 3-10, although highly collimated 
molecular outflows have been detected (NGC 2264G, Lada 
& Fich 1996, Mon R2, Meyers-Rice & Lada 1991, RNO 
43, Padman et al. 1997). Average velocities observed in 
molecular outflows are in the range 10 — 3Gkm s~^, but 
high velocity flows are also observed, with bulk velocities 
up to 100 km s~^ (eg. IRc2, Rodriguez- Franco et al. 1999). 
Also in low- velocity outflows it is possible to find clumps of 
Extremely High Velocity gas (the so-called EHV bullets), 
which shows up as a bump in the tails of spectral lines 
(Cabrit et al. 1997); the mass involved in EHV gas is a 
very small fraction of the total flow mass, from 10"^ to 
10-2 Mq. 

There are a number of properties of the outflows with 
which theoretical models for their acceleration must con- 



front. First we have the distribution of the flow mass 
with velocity, which has typically the form of a power law 
M{v) oc u^, with a break at high velocities. The values 
of 7 range from —1.3 to —2.1 (Cabrit et al. 1997, Mas- 
son & Chernin 1992). The break is found at velocities of 
^25-30 km s^^ and the distribution at higher velocities 
is steeper with slope ~ —3.5, —5. 

We have then the coUimation properties: the flow coUi- 
mation increases with velocity, in fact the higher velocity 
material tends to lie in the center of the outflow, while 
limb-brightening is strongest at low-velocities (L1551- 
IRS5 Uchida et al. 1987, Moriatry-Schieven & Snell 1988; 
NGC 2071 Moriatry-Schieven, Hughes & SneU 1989; NGC 
2024 Richer et al. 1992; NGC 2264G Marguhs et al. 1990). 
In addition molecular outflows show a high degree of bipo- 
larity: in the same lobe the occurrence of overlapping blue- 
shifted and red-shifted emission is rare: the contrast in 
intensity between red and blue-shifted emission increases 
with velocity up to values >20 for high velocity gas (Lada 
& Fich 1996), and this implies a large predominance of lon- 
gitudinal motions (Cabrit et al. 1997, Lada & Fich 1996). 

Many outflows show an apparent linear acceleration, 
i.e. the largest velocities in the molecular material are 
found farthest from the star. Molecular outflows not show- 
ing this property are also observed, for example in NGC 
2024 the velocity is constant over 75% of the outflow (Pad- 
man et al. 1997). 

Finally, Chernin & Masson (1995) studied the distribu- 
tion of mass dM/dz and momentum dP/dz as a function 
of distance along the flow axis and they found that both 
distributions peak in the middle of the lobe, with minima 
near the star and at the extremity of the flow. The same 
trend is observed for the cross sectional area of the flow 
(i.e. the lobes are widest at the location of the momentum 
peak). 

3. Numerical Simulations 

We have studied the entrainment properties of a super- 
sonic jet: the turbulent mixing layer, through which the 
external matter is entrained by the jet motion, is formed 
by the growth and evolution of Kelvin-Helmholtz insta- 
bilities and, therefore, in our simulations we have given to 
the jet an initial small perturbation and we have then fol- 
lowed the flow evolution. The details of the physical and 
numerical setup are given in Paper I, and here we summa- 
rize the main characteristics. Since wc study the evolution 
over long times, we perform the instability analysis using 
the so called temporal approach, in which one studies the 
temporal evolution of the instability in a section of an in- 
finite periodic jet (periodic boundary conditions are used 
at the longitudinal boundaries). The initial velocity and 
density profiles are given by the functional forms 



V^{y,z) = V„ sech 



\nF+z^ 
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Pjot 



t^ — (i^ — 1) sech 



n/? 



where Vg is the velocity on the jet axis, a is the initial jet 
radius, m is a parameter controUing the interface (shear) 
layer width (typically, we set w = 8) and v is the ratio of 
the density at r = yy^ + z^ = cx), p^a, to that on the jet 
axis {y = z ~ 0) a.t t — 0, pjd = Po{y — 0, z = 0). The 
perturbations have the form 



V, 



Vy(x,y,z) = — ^^ sech 



VF+22 



y^^am(nk„x + (f)n), 



(2a) 



y.(x,y,z) = ^sechK^!±i! 
n„ \ a 



cos{nkgX + (t>n) , 



We must notice however that given the nature of our ap- 
proach, i.e. the temporal analysis, we cannot address those 
properties that refer to the distribution of different quan- 
tities along the flow. 

J^.l. Momentum transfer 

The growth of Kelvin-Helmholtz instabilities induces a 
transfer of momentum from the jet to the ambient mate- 
rial, which is then accelerated. The efficiency of this pro- 
cess is very high: Fig. n^ shows the behavior of jet momen- 
tum as a function of time for the three cases of overdense, 
underdense and equal density jet and we see that, in less 
than 20 time units, the jet transfers almost all of its mo- 
mentum to the ambient medium (95% in the case of the 
underdense jet). Our unit of time is the sound crossing 
time over the jet radius that can be expressed as 



tr 



334 ai6 [vsd) yrs 



(26) 



where Vy^^^ = O.OOSV^ is the amplitude of the initial per- 
turbation and {(f>n} are the phase shifts of the various 
Fourier components. The perturbations are thus given by 
a superposition of a number no = 8 of longitudinally pe- 
riodic transverse velocity disturbances, in order to excite 
a wide range of modes. In order to study the entrainment 
and mixing properties we also follow the evolution of a 
tracer, a scalar field, T, passively advected by the fluid. 
In the initial configuration T is put equal to one inside the 
jet and to zero outside, this allows to distinguish the jet 
material from the external material during the whole evo- 
lution. The main control parameters are the Mach number-^ 
of the jet flow, the density contrast between jet and exter- ^ 
nal medium and the actual jet density (for a more detailed 
discussion see Paper I) . The results that we present are for 
a single value of the Mach number M — 10 and three dif- 
ferent values of the density ratio, i.e. for an underdense jet 
(u = 10), for an equal density jet {v = 1) and for an over- 
dense jet {v = 0.1). The conservation equations for mass, 
momentum and energy are integrated using a numerical 
code based on a PPM scheme, on a grid of 256 x 256 x 256 
grid points. The energy equation includes also non equi- 
librium radiative losses, whose detailed form is discussed 
in Paper I. 

4. Results 

In order to determine whether turbulent entrainment that 
develops as a consequence of the evolution of Kelvin-Helm- 
holtz instabilities is a viable mechanism for the accelera- 
tion of molecular outflows, we will now compare the prop- 
erties of the entrained material as they result from our 
simulations with some of the typical properties of molec- 
ular outflows, which have been summarized in Sect. 2. 



where aig is the jet radius a expressed in units of lO^^cm 
and Vse is the sound speed Vg expressed is in units of 10^ cm 
s~^. Therefore, if a jet carries enough momentum flux to 
drive a molecular outflow, this mechanism ensures that 
it would be transferred almost completely to the ambient 
material. 

-Jet rn,omert,tij,Trh 
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Fig. 1. Plot of the jet momentum for heavy (dot-dashed 
line), light (solid line) and cquidense jets (dashed line) 



The next question we can ask is where the momen- 
tum is deposited, more precisely how far from the jet axis 
the ambient material can be accelerated. Due to compu- 
tational limitations, our domain has a limited extension 
in the transversal direction, and therefore we cannot fol- 
low the acceleration of the external material for distances 
larger than 7 jet radii, nonetheless our data allow us to 
give some interesting estimates. We evaluated the rate of 
expansion of the molecular outflow by calculating for each 
time the distance (ij_ from the jet axis within which 50% 
of the deposited momentum is contained (Fig. ||). 
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Fig. 2. Plot of the distance from the jet radius within 
which 50% of ambient medium momentum is found at 
each time, for heavy (dot-dashed line), light (solid line) 
and equidense jets (dashed line) 
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Fig. 3. Logarithmic plot of the mass of ambient medium 
vs its velocity for the light jet case at time t = 18. 
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In the initial times, corresponding to the first stages 
of the instability growth (see Paper I), the ambient ma- 
terial is accelerated only in the immediate surroundings 
of the jet; as time elapses, and the jet starts expanding 
and mixing with the environment, the transversal expan- 
sion of the molecular outflow grows linearly with time, i.e. 
t where the expansion velocity Wexp is similar 



f-^cxp 



for all cases (wcxp — O.lSws for heavy jets, Wcxp = 0.17us 
for equal-density jets, Wcxp = O.lSws for light jets). If we 
extrapolate this linear behavior for longer times we get, 
in a timescale of lO^yrs, a transverse size of the outflow of 
the order of 



d 



ll7. 



4.5 X 10 Us6 ^5 cm 



at the lower end of the observed range; ts is the evolution- 
ary time in units of 10^ yrs. 

4-2. Velocity distributions 

An important constraint that observations provide to the- 
oretical models is the distribution of the mass of the out- 
flow with velocity. The variation of outflow emission (and 
thus of flow mass) is well represented by a power-law shape 
M(y) oc u''' , up to a break velocity, beyond which the mass 
decrease is steeper. 

We computed the distribution of the entrained mass 
with velocity for all our cases, at all times. We found that 
for the light and equal density jet cases, the distribution 
of the ambient mass with velocity can indeed be well rep- 
resented by a power law, with a break at high velocity. As 
an example, we plot, in Fig. 0, the distribution of mass 
with velocity for the light jet case at time t = 18, with the 
two power law fits superimposed. 



7 = 0.09 



-3 r 



-4 r 




0.0 0.5 



1.0 1.5 2.0 2.5 

J'Og(Vj 



3.0 



Fig. 4. Logarithmic plot of the mass of ambient medium 
vs its velocity for the heavy jet case at time i = 6. 

In the heavy jet case, instead, the distribution of mass 
with velocity can be represented by three power-law dis- 
tributions: the mass falls steeply with velocity for low and 
high velocities, while it is almost constant or increases at 
intermediate velocities (Fig. 0). The intermediate flat part 
of the distribution tends to become narrower with time, 
however is still present at the end of our simulation. It 
seems, therefore, that the observed distribution cannot be 
reproduced by the overdense case, at least in the first part 
of its evolution. 

For this reason we concentrated on the analysis of the 
other cases and we examined for them in more detail the 
behavior of the slope of the low velocity region and of the 
position of the break as a function of time, that we plot 
in Fig. H. Two general trends can be recognized: (a) the 
spectral index of the low- velocity region tends to decrease 
with time, i.e. the low-velocity spectrum steepens; (b) the 
break velocity decreases as time elapses. 
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Fig. 5. Trend with time of the index 71 (panel a)) and of 
the break velocity Vb (panel b)) in the power-law distribu- 
tion of ambient mass with velocity m{v) oc v'^; the plots 
refer to the light jet case. 

We can interpret the shape and time evolution of the 
mass distribution in terms of an average velocity profile 
(Stabler 1994), in fact: 



dM 
dv 



dM 
dr 



2'Krp 



dv\ 



(1) 



where M is the mass per unit length and we have as- 
sumed an average homogeneous density distribution. A 
power law velocity profile thus corresponds to a power 
law mass distribution: if v{r) oc r~" we obtain dM /dv oc 
^-1-2/q^ The actual average velocity profiles found in the 
simulations can be well represented by power laws in the 
external low velocity part, while they flatten in the cen- 
tral high velocity region, leading to the steeper part of the 
mass distribution. 

In Fig. g we plot the average velocities for the dense 
jet case at a fixed time, and for the light jet case at two 
selected times. From this figure, taking into account the 
interpretation given above, it is possible to infer the differ- 
ent behaviour in the mass-velocity distribution for heavy 
and light jet cases. For the light jet case (see for example 
the curve corresponding to time i = 18 in Fig. o), the ve- 
locity distribution with radius is fiat at high velocities (0 - 
1.5a) and steeper at intermediate and low velocities (1.5a 
- 3a). In the heavy jet case, instead, we have a flat distri- 
bution at high velocities (0 - 0.5a), followed by a steeper 



distribution at intermediate velocities (0.5a - 2a), while 
at low velocities we have a flattening again, reflecting the 
three components behavior of the mass distribution. 
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Fig. 6. Average distribution of the velocity of the ambient 
medium versus the distance from the jet axis. The solid 
lines are for the light jet case, the dot-dashed line is for 
the heavy jet case. The velocity distribution for the heavy 
jet case has been multiplied by 10 in order to facilitate the 
comparison with the other case. 

Considering the time evolution of the mass distribu- 
tion, we can deduce the temporal variation of the break 
velocity and of the spectral index (Fig. g) from the tem- 
poral variation of the distribution of the average velocity 
with radius. In fact, comparing the two curves in Fig. g 
referring to the light jet case at times i = 16 and i = 18 
(solid lines), we find that the maximum velocity value 
decreases with time, as well as the "break" value where 
the distribution of velocity with radius changes its slope 
{{V) « 2, f = 16; {V) « 1, f = 18). This explains the de- 
crease in the break velocity with time, which is a common 
trend in all the studied cases (Fig. ||, panel b)). Moreover, 
comparing the slope of the average velocity profile at time 
t = 16 with the one at time i = 18 in the low velocity 
range, between 2a and ^ 3a (Fig. 0), we note that the 
velocity distribution at time i = 18 is fiatter since more 
material has been accelerated at larger distances from the 
outfiow axis. This implies a steepening of the mass dis- 
tribution at low velocities, according to Eq. |l|, and this is 
indeed what we see in the light and equal density jet cases 
(Fig. |, panel a)). 

Although the average velocity has the smooth profile 
shown in Fig. O, the actual distribution can however show 
high velocity peaks, which may be correlated with the 
so called "EHV bullets" . We have marked the material 
moving at V^ > 40 km s~^ and integrated over the z di- 
rection. The resulting images are displayed in Fig. M. The 
figure shows that turbulent entrainment is able to accel- 
erate "bullets" of material at high velocities and these are 
found within a few radii from the jet axis. 
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Fig. 7. High velocity structures in the ambient medium, 
integrated over the z direction. 



4-. 3. Collimation 

Observational data on molecular outflows show that the 
flow collimation increases with flow velocity. Moreover, the 
blue-shifted and red-shifted lobes are generally well sepa- 
rated, especially at high velocities. 

To verify whether these features are reproduced in our 
calculations, we considered at each grid point the angle 
between the outflow velocity and the initial jet axis. 
Smaller 9 imply more collimated material. We divided the 
observed velocity range in four intervals, and for each of 
them we calculated the distribution of the ambient mate- 
rial mass with respect to 6. One example of the distribu- 
tions that we obtained is shown in Fig. ||, for the light jet 
case, at time t = 18. 

At low velocity the flow angle is distributed over a wide 
interval, and thus the flow is poorly collimated. The ma- 
terial moving at higher velocities moves at smaller angles 



(a) 



(b) 




S3 0.4 




(c) 



(d) 



0.40 


\ 


0.30 




fv 




J 0.20 






h 


0.10 




\ - 


0.00 




\ 



0.015 


n 


0.010 




J 


- 


0.005 








0.000 


r 


I 



-50 bO 



(a) 0.2 < V,yc, < 1. 2. < V„, < 10. Km s"' 

(b) 1 . < K„/c, < 2. 10.< 7„, < 20. Km s"' 

(c) 2. < K^^/c, < 3.5 20. < V^^t < 35. Km s"' 
(d) V„/c^ > 3.5 V^^, > 35. Km s'' 



Fig. 8. Distribution, in four different velocity intervals, 
of the ambient medium mass with respect to the angle 9 
between the flow velocity and the initial jet axis. These 
plots concern the light jet case, at time t = 18 



respect to the initial flow axis, and is thus more collimated, 
in agreement with observations. This behavior is observed 
in all the studied cases, although a jet moving in a denser 
environment accelerates a more collimated outflow, while 
dense jets accelerate less collimated outflows. This trend 
can be inferred from Table ^ where we report the value 
of the angle 9 within which is contained 90% of the am- 
bient mass whose velocity lies in one of the four intervals 
defined in Fig. y, at three significative times for the three 
cases. The variation with time is in any case quite small. 
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Table 1. Angle containing 90% of the ambient mass mov- 
ing with velocity in the selected ranges, for the three cases, 
at three significative times. 
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Fig. 9. Ratio between the mass of blue-shifted and red- 
shifted material, in the blue-shifted lobe, for outflows mov- 
ing at an angle of 75° (upper panel) and 60° (lower panel) 
respect to the line of view, versus the flow velocity. The 
solid, dashed and dot-dashed lines refer, respectively, to 
the light, equal-density and heavy jet cases. 

We evaluate now how these results compare with the 
observations regarding the overlapping of blue-shifted and 
red-shifted emission in molecular outflows. 

This obviously depends on the angle ip between the line 
of sight and the flow axis (coincident with our x axis) and 
wc choose for our calculation two values for this angle, i.e. 
V' = 75° and ip = 60°. For these two angles we compute 
the ratio of the blue-shifted (positive velocity along the 
line of sight) to the red-shifted (negative velocity along 
the line of sight) mass. In Fig. g we plot this ratio as 
a function of velocity. We notice that the ratio increases 
as this angle diminishes; when the flow is more directed 
towards the observer, in fact, only the material moving 
with large angles respect to the flow axis (i.e. the loosely 
collimated material) will appear red-shifted. 

From our results we notice that the contrast ratio 
between blue-shifted and red-shifted material increases 
with velocity (again, the faster material is the more col- 
limated one), in agreement with the data of, eg. Lada & 
Fich (1996). The contrast ratio, at all velocities, is higher 
for outflows driven by light jets, compared to outflows 
driven by equal density and heavy jets. In Fig. |9| we plot 



the trends for three selected times, and this is possible 
since we did not notice significant differences in the values 
and trends of the contrast ratios at different evolutionary 
times. 

4-.4. Acceleration and dissociation 

One problem of the prompt entrainment mechanism for 
the acceleration of molecular outflows, is that as ambient 
medium is accelerated by the passage of a bow shock, it is 
also dissociated (Downes & Ray 1999), and thus, although 
the bow shock is very efficient in transferring momentum, 
most of it goes to atomic material and the effective fraction 
of momentum transfered by the jet to molecular material 
is low. 
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Fig. 10. Distribution of the ambient accelerated mass as 
a function of its temperature, for different times, and for 
the light jet case. 

In our case we cannot discuss molecular dissociation, 
since we did not include in our computations all the 
physics related to molecular formation and dissociation, 
and to molecular radiative losses. However it is possible 
to have a feeling of what the situation would be, by con- 
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sidering the temperature of the accelerated material. We 
limit the discussion to the light jet case: since we start 
from a configuration where the jet is in pressure equilib- 
rium with the ambient medium, in the :/ = 10 case the 
initial external temperature is 1/10 of the initial jet tem- 
perature (i.e. Text — 1000 K), while in the other cases 
the ambient medium is hotter from the beginning of the 
calculations and cannot be molecular. 

In Fig. M we plot the distribution of the acceler- 
ated ambient material (i.e. external material with velocity 
V > 2 km s~^) with temperature; we note that at inter- 
mediate times (eg. t = 10, 12, 14) a fraction of the ambient 
material is heated, and this is consistent with the fact that 
at this stage strong shocks are driven by the jet into the 
external medium. As shocks become weaker and disap- 
pear, the temperature of the accelerated ambient medium 
decreases. At all times, however, the largest fraction of 
accelerated material maintains a temperature not too dif- 
ferent from the initial value. We can therefore expect that 
the ambient medium accelerated through turbulent en- 
trainment at the jet surface remains mostly molecular. 

5. Summary and Conclusions 

We have shown that turbulent mixing in Kelvin-Helmholtz 
unstable jets is a very efficient mechanism of momentum 
transfer from the jet to the ambient medium: at the final 
stages of the evolution, up to 95% of the jet momentum 
has been transfered to entrained ambient material; since 
the momentum carried by parsec-scale YSO jets is found 
to be comparable to the momentum of typical molecular 
outflows, we infer that this mechanism alone can be in- 
voked to explain the observed outflows. Moreover, oppo- 
site to what happens in bow-shock accelerated material, 
the temperature of the accelerated ambient medium re- 
mains low, and thus momentum is actually acquired by 
molecular material. 

The transverse size of the outflows that are generated 
in this way is located at the lower end of the observed 
range, and our length to width ratio is too high with re- 
spect to the average value for molecular outflows but it 
is similar to that of highly coUimated molecular outflows 
such as, for example NGC 2264G, Mon R2 and RN043. 

An interesting feature of our model is that it repro- 
duces the observed distribution of mass with velocity: it 
has the form of a power law M{v) oc v'* , with a break 
at high velocities. The shape of the distribution and the 
values of the spectral index and of the break velocity are 
particularly well reproduced by the case in which a light 
jet is moving into a denser environment. The most pow- 
erful observed molecular outflows are located in the inner 
regions of molecular clouds, where the ambient medium 
is often so dense that the optical jets are completely ob- 
scured. It is thus plausible that in these cases the atomic 
jets are much lighter than the environment, as implied by 
our model. Conversely, bright optical jets coming out of 



the parent clouds, that would be probably denser than 
the surrounding ambient material, are accompanied by 
very weak molecular outflows, if any (HH 34 and HH 1-2, 
Chernin & Masson 1995). One exception is the bright opti- 
cal jet HH 111, which is associated to a powerful molecular 
outflow (Cernicharo & Reipurth 1996). However, also in 
this case, the blueshifted jet breaking through the surface 
of the cloud is accompanied by a weaker molecular out- 
flow with respect to the redshifted lobe surrounding the 
optically obscured counter-jet (Reipurth & Olberg 1991). 

As far as the distribution of velocity with distance from 
the source is concerned, we could not reproduce the ob- 
served trend: we found that the average velocity of the 
outflow decreases with time (and thus with distance from 
the source, although only qualitative considerations can 
be performed on quantities requiring a translation from 
a temporal to a spatial approach), and this behavior is 
opposite to the observed "Hubble Law" . To explain the 
increase in velocity with distance from the source, many 
possibilities have been invoked: in the frame of the steady 
state model. Stabler (1994) inferred from analytical cal- 
culations that the Hubble law was a consequence of ma- 
terial with progessively higher velocities coming into view 
at greater distances from the star. The Hubble law could 
possibly be due to a decreasing density of the ambient 
medium (Padman et al. 1997), although, more probably, 
high velocities at the head of the flow are actually due to 
the presence of the bow shock that accelerates ambient 
material through prompt entrainment. 

The collimation and bipolarity properties of the out- 
flows generated with Kelvin Helmholtz induced turbulent 
entrainment fit very well the observations: the flow colli- 
mation increases with velocity, and there is a clear sepa- 
ration of the blue-shifted and red-shifted material, which 
increases as well with velocity. An example of a molecular 
outflow for which this behavior has been studied in detail 
is NGC 2264G, Lada & Fich (1996). 

Since these preliminary results are encouraging, future 
work will analyze the combined effects of the two mech- 
anisms, the turbulent entrainment at the jet surface and 
the prompt entrainment at the bow shock, studying the 
propagation of the jet's head in space, and contemporar- 
ily the spatial growth of the Kelvin-Helmholtz unstable 
modes in the jet body. 
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